Artificial gauge fields for the Bose-Hubbard model on a checkerboard 
superlattice and extended Bose-Hubbard model 
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We study the effects of an artificial gauge field on the ground-state phases of the Bose-Hubbard model on a 
checkerboard superlattice in two dimensions, including the superfluid phase and the Mott and alternating Mott 
insulators. First, we discuss the single-particle Hofstadter problem, and show that the presence of a checker- 
board superlattice gives rise to a magnetic flux-independent energy gap in the excitation spectrum. Then, we 
consider the many-particle problem, and derive an analytical mean-field expression for the superfluid-Mott and 
superfluid-altemating-Mott insulator phase transition boundaries. Finally, since the phase diagram of the Bose- 
Hubbard model on a checkerboard superlattice is in many ways similar to that of the extended Bose-Hubbard 
model, we comment on the effects of magnetic field on the latter model, and derive an analytical mean-field 
expression for the superfluid-insulator phase transition boundaries as well. 
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I. INTRODUCTION 

In the usual Bose-Hubbard model f^, competition between 
the kinetic and potential energy terms leads to two phases: a 
Mott insulator (MI) when the kinetic energy is much smaller 
than the potential energy and a superfluid (SF) otherwise. The 
incompressible MI phase has an excitation gap so that the in- 
coherent bosons are localized, and that a slight change in the 
chemical potential does not change the number of bosons on 
a particular lattice site. The compressible SF phase, however, 
is gapless, and the coherent bosons are delocalized over the 
entire lattice. 

Recent advances with ultracold bosonic atoms loaded into 
optical lattices have made it possible to simulate Bose- 
Hubbard type many-particle Hamiltonians in a tunable setting. 
For instance, the ability to control on-site boson-boson inter- 
actions has paved the way for observing SF and MI phases as 
well as the transition between the two iH 0] • In addition, a 
new technique has recently been developed that allowed the 
production of fictitious magnetic fields which can couple to 
neutral bosonic atoms IH U. These fictitious magnetic fields 
are produced through an all optical Raman process, couple to 
a fictitious charge, but produce real effects like the creation of 
vortices in the SF state of bosons. Such an ability to control 
the strength of the fictitious magnetic fields combined with 
the ability to control the strength of the interparticle interac- 
tions may allow exploration of new phenomena in the near 
future lHHil. 

In contrast to its simplicity, the Bose-Hubbard model is not 
exactly solvable even in one dimension. Therefore, it is de- 
sirable to have a much simpler toy-model which exhibits all 
the salient properties of the Bose-Hubbard model, while also 
being more amenable to analytical treatment. One of the most 
prominent candidates is the hardcore Bose-Hubbard model on 
a checkerboard superlattice, for which the existence of SF and 
MI phases at half filling in three dimensions |13], as well 
as a direct transition between the two IJ4iil have rigorously 
been shown. In addition, this model and its correlation func- 
tions are exactly solvable in one dimension [Ts], thanks to 
the existence of a mapping between the hardcore bosons and 



noninteracting fermions. In the absence of a magnetic field, 
we have recently analyzed the ground-state phase diagram 
of this model in one, two and three dimensions using mean- 
field approximation and strong-coupling expansion, and com- 
pared them with the numerically exact results obtained from 
the stochastic series expansion algorithm followed by finite- 
size scaling 03- Giveijjhat checkerboard superlattices have 
already been realized [19] using multiple wavelength laser 
beams, we extend previous works in two important directions. 
First, we relax the hardcore constraint and study the ground- 
state phase diagram of the softcore Bose-Hubbard model on a 
checkerboard superlattice, via mean-field decoupling approx- 
imation. Second, we study the effects of uniform magnetic 
field on the ground-state phase diagram in two dimensions. 

The rest of this paper is organized as follows. First, we re- 
view the model at hand in Sec. [Ill and present a qualitative 
description of its phase diagram. Then, we study the single- 
particle Hofstadter problem in Sec.|IIIl and show that the pres- 
ence of a checkerboard superlattice gives rise to a magnetic 
flux-independent energy gap in the excitation spectrum. The 
many-particle problem is discussed in Sec. lIVI where we de- 
rive analytical expressions for the SF-insulator phase transi- 
tion boundaries within the mean-field decoupling approxima- 
tion. Then, we give a brief discussion and summary of our 
results in Sec. |Vl Finally, we conclude the paper with Ap- 
pendix A, where an analytical mean-field expression for the 
SF-insulator phase transition boundaries is derived for the ex- 
tended Bose-Hubbard model. 



II. HAMILTONIAN 



In this paper we study the effects of magnetic field on the 
Bose-Hubbard model on a checkerboard superlattice. For this 
purpose, we consider a two-dimensional square lattice de- 
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scribed by the Hamiltonian 

H = - ^%e*^-ajaj + ^^Mri^- 1) 

ij i 

-C^(-l)'^-n,-M^H„ (1) 

i i 

where a\ (ui) creates (annihilates) a boson on site i and 
rii = ajfli is the on-site boson number operator The hopping 
matrix is assumed to connect two nearest-neighbor lattice 
sites (tij = t for i and j nearest neighbors and otherwise) 
belonging to different sublattices, e.g. the even sublattice A 
and the odd one B,U > Ois the strength of the on-site boson- 
boson repulsion, ^ is the chemical potential, and C > is the 
amplitude of the alternating checkerboard superlattice poten- 
tial such that (Ti = (1) on sublattice A (B). The phase factor 

9ij = (l/0o) // -A-o(r) • dr takes into account the effects of 
a uniform magnetic field that is applied perpendicular to the 
lattice, where Ao(r) is the vector potential and (po = hc/e 
is the magnetic flux quantum. All of our results recover the 
nonmagnetic ones when the magnetic flux tends to 0. 

Let us first analyze the atomic {t ~ 0) limit of this Hamil- 
tonian at zero-temperature. In this limit, there is no kinetic 
term, and the boson number operator commutes with the 
Hamiltonian, so every lattice site is occupied by a fixed num- 
ber rii = (rii) of bosons. Here, (...) is the thermal average, 
and the average boson occupancy rii is chosen to minimize 
the ground-state energy for a given /x. In particular, when 
C = 0, this model is translationally invariant, and the ground- 
state boson occupancy is the same for all sites. 

For instance, in the hardcore boson (U oo) limit lfT6ll . 
while the lattice is completely empty for /i < and the min- 
imal energy configuration corresponds to a vacuum of parti- 
cles or a hole band insulator (VP, since = for all i), it 
is completely full for fj, > and the minimal energy config- 
uration corresponds to a vacuum of holes or a particle band 
insulator (VH, since = 1 for all i). The ground-state en- 
ergy of these phases is degenerate at /i = 0. However, when 
C > 0, the ground state has an additional half-filled insulating 
(incompressible) phase characterized by a crystalline order in 
the form of staggered boson densities. For the nearest neigh- 
bor lattice sites i and j, (rii) = ua ^ I for the sublattice A 
and (rij) — ns = for the sublattice B. This phase resides in 
the region < C, and it is sandwiched between the VP and 
VH phases. Since the checkerboard superlattice breaks the 
translational invariance of the lattice, it directly causes such 
an alternating density pattern. For this reason, this phase is 
often called a MI 1 14, 16], to distinguish it from a true charge- 
density-wave (CDW) phase, for which the translational invari- 
ance is broken spontaneously due for instance to the presence 
of nearest-neighbor interactions. (See Appendix A). In this 
paper, for simplicity we call this alternating density pattern an 
alternating MI (AMI) with (1, 0) fillings to prevent confusion. 

On one hand, for softcore bosons with U > 2C ^ 0, the 
ground state alternates between the AMI and MI phases as a 
function of increasing ^, where the chemical potential widths 
of AMI and MI lobes ai-e 2C and U - 2C, respectively (see 
Fig.|4|. On the other hand, for softcore bosons with 2C > U, 



the ground state has only AMI insulators. For instance, when 
2U > 2C > U, the ground state is a VP with (0,0) fiflings 
for /I < — C; it is an AMI with (1,0) fillings for /i between 
-C and U ~C;\i is an AMI with (2, 0) fiflings for ^ between 
U - C and C; it is an AMI with (2, 1) fiflings for ^ between 
C and 2U - C; it is an AMI with (3, 1) fiflings for /i between 
2U — C and U + C, and so on. As t increases, the range of /i 
about which the ground state is insulating decreases, and the 
MI and AMI phases disappear at a critical value of t, beyond 
which the system becomes compressible. We note that, unlike 
the compressible SF phase of the usual Bose-Hubbard model, 
the compressible phase in this model is more like a supersolid 
(SS), where the SF and AMI orders coexist even for arbitrar- 
ily small C. However, since the checkerboard superlattice 
breaks the translational invariance of the lattice, we call the 
compressible phase of this model a SF to distinguish it from a 
true SS for which the translational invariance is broken spon- 
taneously. Having discussed the atomic limit, let us now dis- 
cuss the noninteracting single-particle energy spectrum when 
t ^ 0. 

III. SINGLE-PARTICLE PROBLEM 

In the absence of a checkerboard superlattice, i.e. when 
C = 0, the single- part icle excitation spectrum is the usual 
Hofstadter butterfly ||20|] . and here we generalize it to C 7^ 0. 
For this purpose, we choose the Landau gauge for the vector 
potential [Ao(r) = (0, Bqx, 0)], which leads to a uniform Bq 
field in the z direction, and the strength of the magnetic field 
Bf) is related to the magnetic flux $ via $ = B^l^. Denoting 
the coordinates of lattice sites by i = (x = nl^y — ml), 
this gauge simply implies 9ij = for hoppings along the x 
direction, i.e. between {n,Tn) and {n ± l,m); and 9ij = 
±2TT(t)n for links along the y direction, i.e. between (n, m) 
and (71, m ± 1), where (f) — $/(27r(/)o). In the Landau gauge, 
the Hamiltonian given in Eq. ([T]) can be written as 

nm 

-C^(-l)"+™aLa«™, (2) 

nm 

where 77 + m is even (odd) for sublattice A (B). Here, we set 
U = and n — for the single-particle problem. 

In the Landau gauge, taking (j> — p/q, where p and q are 
integers with no common factor, while the Hamiltonian given 
in Eq. (|2]i maintains its checkerboard translational invariance 
in the y direction, i.e. it remains the same under 2 steps 
(m 171 + 2), it requires q* steps for translational invariance 
in the x direction. For the square lattice considered, the period 
q* = q (2g) for even (odd) q values. Therefore, the Bloch 
theorem tells us that the 1 st magnetic Brillouin zone is deter- 
mined by -7r/2 < ky£ < Tr/2 and —n/q* < k^i < Tr/q*. 
This increased periodicity motivates us to work with a super- 
cell of 2 X g* sites (in x and y directions, respectively) as 
shown in Fig.[Tl 

The single-particle excitation spectrum is de- 
termined by solving the Schrodinger equation 
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FIG. 1. (Color online) The 2 x q* supercell is shown where q* — q 
{2q) for even (odd) q values. 



Hsc'i'sc — Ec{(f>)'^sc for all k values in the 1st 
magnetic Brillouin zone. Choosing the wavefunction, 

vl',, = ii;^ , V-f , i^i , <. , ^f, where 

ijjn^'^^ denotes the nth site of sublattice A (B), the 2q* x 2q* 
matrix 
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(3) 



describes the supercell with periodic boundary (Bloch) con- 
ditions. Here, a„ — —2tcoa{kyi + 2TTnp/q). Equation (|3]l 
is the generalization of the Hofstadter problem to the case of 
a checkerboard superlattice, and it reduces to the usual result 
whenC = |20]. 

In Fig. 121 we show the single-particle excitation spectrum 
Ec{4>) as a function of the magnetic flux (j) = p/q for two 
values of C: (a) C = and (b) C = i. In both figures, Ec{<j}) 
is shown to be symmetric around p/q = 1/2, i.e. Ec{(t>) — 
Ec{^ — (f>)- This is simply because the magnetic flux values 
that add up to 2TT(j)o are equivalent, e.g. = and 0=1. This 
also means that the maximal magnetic field Bq that can be 
applied corresponds to p/q = 1/2. A list of minimal single- 
particle excitation energies e(0) = mint Eo{(t)) jt of the usual 
Hofstadter butterfly is also given for a number of pjq values 
in Table.!] 

In contrast to the presence of zero-energy excitations for 
all possible values in Fig. |2ja), the most important differ- 
ence in Fig.|2b) is the presence of a flux-independent energy 
gap. Our numerical calculations show that the energy gap is 
exactly 2C. Since the on-site energy difference between the 
sublattice A and B is 2C, this result is not so surprising at least 
in the ^> limit. In an earlier work |16], in the absence 
of a magnetic field, we showed that the single-particle exci- 
tation spectrum becomes £'c(0) = i-^/C^ + -Eq (0), where 
-Eo(O) = — 2tcos(fcx^) — 2tcos{ky£) is the usual single- 
particle excitation spectrum. Similar to the nonmagnetic case, 
our numerical calculations show that 



between the kinetic and potential energy terms of the many- 
particle Hamiltonian when t ^ 0. 



Ec{<P) ^ ±\ + Eli^) 



(4) 



holds exactly in the presence of a magnetic field, where 
Eq {(p) — Ec=o (</>) is the usual single-particle excitation spec- 
trum of the Hofstadter butterfly. Having discussed the single- 
particle problem, now we are ready to analyze the competition 



TABLE I. A list of minimal single-particle excitation energies 
e{<j)) = mink-Eo(0)/t of the usual Hofstadter butterfly (C = 0) is 
given for a number of magnetic flux (f> — ^/[2-n(j)o) — p/q values. 



p/q 


^(0) 


p/q 


^(0) 


p/q 


^(0) 


p/q 


^(0) 


1/2 


-2.828 


2/3 


-2.732 


3/4 


-2.828 


4/5 


-2.966 


1/3 


-2.732 


2/5 


-2.618 


3/5 


-2.618 


4/7 


-2.611 


1/4 


-2.828 


2/7 


-2.725 


3/7 


-2.611 


4/9 


-2.630 


1/5 


-2.966 


2/9 


-2.881 
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4/11 
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1/6 


-3.096 


2/11 


-3.028 


3/10 


-2.698 


4/13 
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1/7 
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2/13 
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3/11 
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4/15 
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2/15 
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3/13 
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FIG. 2. (Color online) The single-particle excitation spectrum 
Ec {4>) of the Hamiltonian given in Eq. ([3} (in units of t) is shown as 
a function of the magnetic flux — <^/{2-n(j)o) — p/q for two values 
of C: (a) C = and (h) C = t. Note that a flux-independent energy 
gap opens in (b), which exactly equals to 2C. 



IV. MANY-PARTICLE PROBLEM 



For the many -particle problem, we add the boson-boson in- 
teraction and chemical potential terms to Eq. dU, and obtain 



nrn 

+ / ^ O'nm^nmjO'nm^nrn ^ Ij 
nm 

^C^^(— 1) QnmQ"™ ~ ^lim^nm- (5) 

nm nm 

For illustrative purposes, let us first study the ground-state 
phase diagram of this model for the hardcore bosons, for 
which the calculation becomes considerably simpler com- 
pared to the softcore case, but yet nontrivial. 



A. Hardcore bosons 

The hardcore boson operators satisfy the constraint a^^^^ — 
'^nm — 0' which prohibit multiple occupancy of lattice sites, 
as dictated by the infinitely large on-site boson-boson repul- 
sion (?7 — > oo). In this limit, the many-particle Hamiltonian 
given in Eq. Q becomes invariant under the transformation 
flnm alii.m or a,\ „j^;^, which corresponds to a 

shift of one lattice site in a; or y direction. This symmetry op- 
eration, which can be immediately read off from the Hamilto- 
nian, corresponds to a particle-hole exchange combined with 
swapping A and B sublattices. It leads to a — > — /i symme- 
try, i.e. the phase diagram is symmetric around /i = 0. 

It turns out that the exact SF-VP and SF-VH phase tran- 
sition boundaries can be easily obtained analytically. The 
simplest argument leading to this conclusion stems from the 
fact that this boundary is determined by the addition of a sin- 
gle particle (hole) to the completely-empty (-filled) lattice. It 
can then be argued that whether one is dealing with hardcore 
bosons or noninteracting spinless fermions makes no differ- 
ence, as the particle statistics plays no role. This further means 
that one needs only to diagonalize the single-particle Hamil- 
tonian and find the energy difference between the completely- 
empty (-filled) lattice and the state with one particle (hole). 
The single-particle spectrum is already given in Eq. (|4|, and 
this procedure leads to 

/i - ±v/C2+ £2(0)^2, (6) 

where e{4)) = minki?o(0)/^ are the minimal single-particle 
excitation energies of the usual Hofstadter butterfly which are 
shown as big blue dots in Fig.|2a). In other words, the minus 
(plus) branch in Eq. ^ is determined by the minimal (maxi- 
mal) single-particle excitation energies of the Hofsdtadter but- 
terfly when C 7^ 0. The minimal single-particle excitation en- 
ergies are shown as big blue dots in Fig.|2b). We emphasize 
that Eq. (|6|l is exact for two-dimensional square lattices, and 
the minus (plus) sign determines the SF-VP (-VH) phase tran- 
sition boundary. Note again that Eq. ^ reduces to the known 
result for the nonmagnetic case II61I in the e(0 — > 0) = —4 
limit. 

Unhke the SF-VP and SF-VH phase boundaries, the SF- 
AMI phase transition boundary cannot be determined exactly, 
since the exact many-particle wave function for the AMI 
state is not known. However, this can be achieved via the 
mean-field decoupling approximation. Within this approxi- 
mation, the hopping terms in the Hamiltonian given in Eq. (|5]l 
are decoupled according to, al^an+i,„i = (at^)a„+i^,„ + 
al^{an+i,m) - (at„)(a„+i,„), where the expectation val- 
ues (pnm{4>) — {O'nm) Correspond to the mean-field SF order 
parameters. Note that there are at most 2q*, i.e. total number 
of sites in a supercell, distinct ^pnm and the SF (AMI) phase 
is determined by the nonzero (zero) value of any one (all) of 
them. In particular, there are two distinct SF order param- 
eters (one for each sublattice) in the nonmagnetic case even 
for arbitrarily small C, as long as C ^ 0. We checked that 
these order parameters differ from each other for all parame- 
ter space via Gutzwiller ansatz calculations. Therefore, unlike 
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the compressible SF phase of the usual Bose-Hubbard model, 
the SF phase in the checkerboard model is more like a SS, 
where the SF and AMI orders coexist. 

Performing a second-order perturbation theory in ipnrn 
around the MI phase, and following the usual Landau pro- 
cedure for second-order phase transitions, i.e. minimizing the 
ground-state energy as a function of (pnm, we eventually ar- 
rive at the phase transition boundary equation 



(7) 



where the plus (minus) sign corresponds to particle (hole) ex- 
citations above the AMI phase. We note that although ipnm 
are gauge dependent, the phase transition boundary itself is 
not 10,13,1^. Alternatively, Eq. ^ follows directly from 
the strong-coupling expansion of the ground-state energy of 
the AMI phase with respect to the hopping term 1 16], where 
the first nontrivial hopping dependence of the phase transition 
boundary arises from the maximal eigenvalue of the T = t • t 
matrix 12V\. Here, the elements of = 'J2k^'-ktkj are 
such that T^jfj = e^(0)i^/i- Note again that Eg^^© re- 
duces to the known result for the nonmagnetic case 11611 in the 
£((/) ^ 0) = -4 limit. 




II /C 

FIG. 3. (Color online) The ground-state phase diagram of hardcore 
bosons is shown as a function of the chemical potential jj, (in units of 
C), hopping 4:t (in units of C) and magnetic flux (j) — <^/{2ncj>o) — 
p/q, as obtained from Eqs. ^ and Q. The phase diagram is sym- 
metric around fi = and p/q — 1/2 as explained in the text, and the 
intriguing structure of the phase boundaries (which is clearly seen 
in the contour plot) is due to the minimal single-particle excitation 
energies e{<j}) shown as big blue dots in Fig.[2{a). 

In Fig. [3] we show the ground-state phase diagram as a 
function of the chemical potential fi, hopping 4t and mag- 
netic flux (f) = p/q, that is obtained from Eqs. (|6]l and (|7]). 
As we argued above, the phase diagram is symmetric around 
/i = and p/q = 1/2. The latter symmetry is in agreement 
with the earlier findings on the usual Bose-Hubbard model 
01 ■ What is more interesting is the intriguing structure of the 
phase transition boundaries on the minimal single-particle ex- 



citation energies e ((/)). In addition, the incompressible (com- 
pressible) AMI (SF) phase grows (shrinks) when the magnetic 
field increases from zero, due to the localizing effects of the 
magnetic field on bosons. All of these observations are similar 
to earlier findings on the usual Bose-Hubbard model IH-Ql. 



B. Softcore bosons 

Having studied the effects of magnetic field on the hardcore 
bosons, let us now analyze the ground-state phase diagram 
of the Hamiltonian given in Eq. Q for the softcore bosons. 
Since the exact many-particle wave functions for the AMI and 
MI states are not known, we again obtain the phase diagram 
via the mean-field decoupling approximation. Following the 
recipe given in the previous section, the SF-MI and SF-AMI 
phase transition boundaries are found to be determined by 
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HA + l 



riB 



1 



riB 



Uuf 



(8) 



C-Ai UiriB - I) + C - 

which gives a quartic equation for /i. This equation is vaUd 
for all C, and e{(p) — mint i?o(</')/i are the minimal single- 
particle excitation energies of the usual Hofstadter butterfly. 
In the hardcore (U — s> oo) limit, note that Eq. (O reduces 
to Eq. O when ua — 1 and = 0, and to Eq. (|6]l when 
nA = ns = or riyi — ns = 1. In addition, Eq. (O 
reduces to the known result iTfl for the usual Bose-Hubbard 
model when ua = ns = no and C = 0, and it also agrees 
with the recent numerical calculations iflTiflsIl in the absence 
of a magnetic field. Since a simple closed form analytic solu- 
tion for /i is not possible when C ^ 0, we solve Eq. (O with 
MATHEMATICA for each of the AMI and MI lobes sepa- 
rately. 

In Fig.lH we set C = 0.2U and show the ground-state phase 
diagram as a function of the chemical potential /i and hop- 
ping At for p/q — 1/1 (equivalent to zero magnetic field) and 
for p/q = 1/2 (maximum magnetic field). As discussed in 
Sec. |II] the ground state alternates between the AMI and MI 
phases as a function of increasing /i. The chemical potential 
widths of AMI and MI lobes are OAU and 0.6U, respectively, 
but the size of the AMI (MI) lobes increase (decrease) as a 
function of increasing C/U (not shown), since a nonzero C is 
what allowed AMI states to form in the first place. In addition, 
the incompressible (compressible) AMI and MI (SF) phases 
grow (shrinks) when the magnetic field increases from zero, 
due again to the localizing effects of the magnetic field on 
bosons. Since the phase diagram of the Bose-Hubbard model 
on a checkerboard superlattice is in many ways similar to that 
of the extended Bose-Hubbard model, we discussed the latter 
model in Appendix A. 



V. CONCLUSIONS 

In this paper, we studied ground-state phases of the Bose- 
Hubbard model on a checkerboard superlattice in two dimen- 
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FIG. 4. (Color online) The ground-state phase diagram of softcore 
bosons is shown as a function of the chemical potential \i (in units of 
V) and hopping 4t (in units of C/), as obtained from Eq. {Sjl. Here, 
we set C = 0.2(7, and pjq = 1/1 (1/2) is shown as solid red 
(dashed black), corresponding to zero (maximum) magnetic field. 
The AMI and MI (SF) phases grow (shrinks) when the magnetic field 
increases from zero, due to the localizing effects of the magnetic field 
on bosons. 



sions. First, we discussed the single-particle Hofstadter prob- 
lem, and showed that the presence of a checkerboard super- 
lattice gives rise to a magnetic flux-independent energy gap 
in the excitation spectrum. Then, we considered the many- 
particle problem, and derived analytical mean-field expres- 
sions for the SF-MI and SF-AMI phase transition boundaries. 
We showed that the size of incompressible insulator phases 
grow when the magnetic field increases from zero, due to the 
localizing effect of the magnetic field on bosons. In addi- 
tion, since the phase boundaries are functions of the minimal 
single-particle excitation energies, they have an intriguing de- 
pendence on the magnetic flux. 

We also showed that the phase diagram of the Bose- 
Hubbard model on a checkerboard superlattice is in many 
ways similar to that of the extended Bose-Hubbard model. 
In particular, the compressible phase in the former model is 
more like a SS, where the SF and AMI orders coexist even 
for arbitrarily small C. However, since the checkerboard su- 
perlattice breaks the translational invariance of the lattice, we 
call the compressible phase of this model a SF to distinguish it 
from a true SS for which the translational invariance is broken 
spontaneously. For completeness, we discussed the effects of 
magnetic field on the the extended Bose-Hubbard model as 
well, and derived an analytical mean-field expression for the 
SF-MI and SF-CDW phase transition boundaries. 



In this paper, we relied on the mean-field theory which 
is known to be sufficient in describing only the qualitative 
features of the phase diagram when there is no magnetic 
field. In particular, for the hardcore Bose-Hubbard model on 
a checkerboard superlattice, we have recently shown that the 
mean-field theory has a large quantitative discrepancy from 
the numerically exact Quantum Monte Carlo results especially 
in lower dimensions iflal . However, due to the infamous sign 
problem, such numerical calculations cannot be performed in 
the presence of a magnetic field, and therefore, it is not clear 
to us whether the mean-field theory is sufficient in this case or 
not. We hope that the intriguing structure of the phase transi- 
tion boundaries on the minimal single-particle excitation en- 
ergies of the Hofstadter butterfly, pedicted by the mean-field 
theory, could be observed in the experiments or verified via 
other exact means in the future. 

We plan to extend this work at least in one important di- 
rection. There is some evidence that the ground-state phase 
diagram of the extended Bose-Hubbard model includes a SS 
phase, in dimensions higher than one I22l- l24l l26ll . The lo- 
calizing effect of magnetic field on such a phase is yet to be 
studied, and it is not obvious whether SS region would grow or 
shrink when the magnetic field increases from zero. Although 
numerical calculations based on the Gutzwiller ansatz are ex- 
pected to give gauge dependent results for the SS-SF phase 
boundary because of the mean-field nature of the ansatz, they 
would provide a good qualitative insight into this problem. 
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Appendix A: Extended Bose-Hubbard model 

In many ways, the phase diagram of the Bose-Hubbard 
model on a checkerboard superlattice turned out to be simi- 
lar to that of the extended Bose-Hubbard model. Therefore, in 
this appendix, we comment on the effects of uniform magnetic 
field on the insulating phases of the latter model. In contrast to 
our model where the translational invariance is broken due to 
checkerboard superlattice, the translational invariance is bro- 
ken spontaneously in the extended model, leading to CDW 
modulations. The extended Bose-Hubbard Hamiltonian with 
the on-site ([/ > 0) and nearest-neighbor (V > 0) boson- 
boson repulsions can be written as 



(Al) 
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Here, we again consider a two-dimensional square optical lat- 
tice. For U > AV ^ 0, it is well-known that the ground state 
has two types of insulating phases Il22l - l25ll . The first one is the 
MI phase where, similar to the usual Bose-Hubbard model, 
the ground-state boson occupancy is the same for every lat- 
tice site, i.e. (n^) — uq. The second one is the CDW phase 
which has crystalline order in the form of staggered boson oc- 
cupancies, i.e. (rii) = ua and (n^ ) = ub for i and j nearest 
neighbors. 
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FIG. 5. (Color online) The ground-state phase diagram of the ex- 
tended Bose-Hubbard model is shown as a function of the chemical 
potential ^ (in units of U) and hopping At (in units of U), as obtained 
from Eq. Here, we set V = 0.1(7, and p/q = 1/1 (1/2) is 

shown as solid red (dashed black), corresponding to zero (maximum) 
magnetic field. 

When U > AV, in the atomic t — {) limit, the ground state 
alternates between the CDW and MI phases as a function of 
increasing /i, where the chemical potential widths of CDW 



and MI lobes are AV and U, respectively. As t increases, 
the range of fi about which the ground state is insulating de- 
creases, and the MI and CDW phases disappear at a critical 
value of t, beyond which the system becomes compressible 
(SF or SS). On the other hand, the ground state has only CDW 
insulators when U < AV. The chemical potential width of all 
CDW insulators is U, and the ground state is a CDW insu- 
lator with (riOjO) fillings for {uq — 1)U < /i < noU. As 
t increases, the CDW phases disappear at a critical value of 
t, beyond which the system first becomes a SS then a SF at 
a much larger t with a very large region of SS phase ll26ll . 
Within the mean-field decoupling theory, the phase transition 
boundaries are determined by 



UA + l 



riA 



UuA + AVnB-H U{nA - 1) + AVub ~ fj, 
Ub + 1 Ub 



UnB+AVuA-l-i U{nB - I) + AVuA- fJ. 

(A2) 



which gives a quartic equation for /i. Here, e((/)) — 
mink Eq {(f)) /t depends on the minimal single-particle exci- 
tation energies of the usual Hofstadter butterfly. Note that 
Eq. (IA2l i reduces to the known result for the nonmagnetic 
case 125.1 in the e(0 — >■ 0) = —4 limit, and it reduces to 
the known magnetic result lUl-Ql when ha — tlb — tlq and 
V = 0. Equation iA2i shows that the MI lobes are sepa- 
rated by AV, but their shapes are independent of V within the 
mean-field decoupling approximation; in particular, the criti- 
cal points for the MI lobes are independent of V. 

In Fig.|5] we set y = 0. It/ and show the ground-state phase 
diagram as a function of the chemical potential /i and hop- 
ping At for p/q =1/1 (no magnetic flux) and for p/q =1/2 
(maximum magnetic flux). As discussed above, the ground 
state alternates between the CDW and MI phases as a func- 
tion of increasing fi. The chemical potential widths of CDW 
and MI lobes are O.AU and U, respectively, but the size of 
the CDW (MI) lobes increase (decrease) as a function of in- 
creasing V/U (not shown), since a nonzero V is what allowed 
CDW states to form in the first place. In addition, the in- 
compressible (compressible) CDW and MI (SF or SS) phases 
grow (shrinks) when the magnetic field increases from zero, 
due to the localizing effects of the magnetic field on bosons. 
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